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Abstract 



We describe the dynamics of fiuxons moving in a frustrated Josephson 
junction with p, d, and /-wave symmetry and calculate the I — V character- 
istics. The behavior of fiuxons is quite distinct in the long and short length 
junction limit. For long junctions the intrinsic flux is bound at the center 
and the moving integer fluxon or antifluxon interacts with it only when it 
approaches the junction's center. For small junctions the intrinsic flux can 
move as a bunched type fluxon introducing additional steps in the I — V 
characteristics. Possible realization in quantum computation is presented. 
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I. INTRODUCTION 



The determination of the order parameter symmetry in high-T c superconductors is a 
problem which has not yet been completely solved [ffl-E). The Josephson effect provides a 
phase sensitive mechanism to study the pairing symmetry of unconventional superconduc- 
tors. In Josephson junctions involving unconventional superconductors, the sign change of 
the order parameter with angle measured from the x-axis in the ab plane introduces an 
intrinsic phase shift of ir in the Josephson current phase relation or alternatively a negative 
Josephson critical current. The effect of shifting the phase by ir is equivalent to the shifting 
the critical current versus the magnetic flux pattern in a squid that contains a 7r junction 
(called frustrated junction) by $o/2, where <3>o = h/2e is the flux quantum 0. 

The presence of spontaneous or trapped flux is a general property of systems where a 
sign change of the pair potential occurs in orthogonal directions in fc-space. Its existence has 
been predicted for example in ruthenates [0] where the pairing state is triplet as indicated 
in the Knight shift measurements || and the time-reversal symmetry is broken as shown by 
the muon spin rotation (fiSR) experiment where the evolution of the polarization of the 
implanted muon in the local magnetic environment of the superconductor gives information 
about the presence of spontaneous magnetic field. Moreover the pairing state has line nodes 
within the gap as indicated by the specific heat measurements [ID|. This spontaneous flux 
shows a characteristic modulation with the misorientation angle within the RuCVplane that 
can be checked by experiment [|Il[ . 

The one dimensional Josephson junction with total reflection at the end boundaries, 



between s-wave superconductors, supports modes of resonant propagation of fluxons [12 



In the plot of the current- voltage (I — V) characteristics these modes appear as near-constant 
voltage branches known as zero field steps (ZFS) Hi3HT5|j. They occur in the absence of any 
external field. The ZFS appear at integer multiple of V\ = $oCs/Z, where cs is the velocity of 
the electromagnetic waves in the junction, and / is the junction length. The moving fluxon 
is accompanied by a voltage pulse which can be detected at the junction's edges. 
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When the contact between a and it junctions, that contains an intrinsic half-fTuxon, 
is current-biased the half-fluxon becomes unstable for certain values of the external current 



with respect to transforming into an anti-half-fluxon and emitting an integer ffuxon 16 



Also when a — n — junction, that contains two half vortices, is current-biased, for certain 
critical current, a transition occurs between the two degenerate ffuxon configurations and a 
voltage pulse is generated [[[7} . 

In this paper we study the dynamic properties of fluxons and calculate the I — V charac- 
teristics in frustrated junctions with Bi g , E u , B\ g x E u pairing symmetry. The last two are 
candidates pairing states for ruthenates f7j. The nodeless p-wave order parameter with E u 
symmetry has been proposed by Rice and Sigrist |1| while the Bi g x E u has been proposed 
by Hasegawa et al. ||19|| . In junctions involving unconventional superconductors the behavior 
of fluxons is typically different in the long and short length junction limit. In the long limit 
the fractional fluxon is confined at the center and the moving fluxon interacts with it only 
when it approaches the center. However in the short limit the bound fluxon becomes able 
to move as a bunched type solution with integer or half integer magnetic flux. For the Bi g 
case the I — V pattern is shifted by a voltage that corresponds to the intrinsic phase shift. 
Also the frustrated Josephson junction can be considered as a way to build a quantum 'bit' 
(qubit) which is the generalization of the 'bit' of the classical computer. 

The article is organized as follows. In Sec. II we develop the model and discuss the 
formalism. In Sec. Ill we present the results for the long junction and in Sec. IV for the 
shorter junction limit. In Sec. V we discuss the implementation of the qubit and finish with 
the conclusions. 



II. CORNER JUNCTION MODEL 

We consider the junction shown in Fig. 1(a) between a superconductor A with a two com- 
ponent order parameter and a superconductor B with s-wave symmetry. The supercurrent 
density can be written as: 
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J{4>) = J c sin(0 + c ), (1) 

where J c is the Josephson critical current density, is the relative phase difference between 
the two superconductors and 4> c is the intrinsic phase shift. We describe a frustrated junction 
of length I, i.e., the two segments have different characteristic phases, i.e., <\> c \ in < x < | 
and C 2 in | < x < I. By introducing an extra relative phase in one part of this junction, 
this one dimensional junction can be mapped in the corner junction that is seen in Fig. [IJ 
(b). The characteristic phases <p c i and C 2 distinguish the various pairing symmetries and 
can be seen in table |. For the orientation of the junction that we consider in which the 
a and b crystal axes are at right angles to the interface a simple calculation llfllllf gives 
J c = 1 for the pairing states that we consider. 

The phase difference across the junction is then the solution of the time dependent sine- 
Gordon equation 

d 2 (f> d 2 (j) d(j) 
with the following inline boundary conditions 

Tx lx=0 ' 1 = 4' (3) 



where the time t is in units uJq 1 , where ujq = is the Josephson plasma frequency, C 

is the capacitance per unit length. 7 = G/ujqC is the damping constant which depends on 
the temperature, G _1 is an effective normal resistance. The value used in the numerical 
calculations is 7 = 0.01. The length x is normalized in units of the Josephson penetration 



depth A j = ./— 2s — , L p is the inductance per length and is given by L p = fi d, where 

\l 27T ,JcLp 

d = 2Xl + t ox is the magnetic thickness of the junction, is the London penetration depth, 
t ox is the thickness of the insulating oxide layer, and fiQ = 4tt10~ 7 H/m. The velocity of the 



electromagnetic waves in the junction is given by c$ = ^L p C. I is the normalized inline 
bias current in units of Xj Jo- 
in previous publication we used overlap boundary conditions, where the current is 
uniformly distributed in space. However in actual experiments in the s-wave case, the 



biased current in the overlap geometry may be concentrated at the edges within a length 



A j rather than distributed in space J2U|. Therefore it is more appropriate to use inline 
boundary conditions. However in the case of the overlap geometry only details of the ffuxon 
propagation are quantitatively different e.g. the oscillations of the bound fluxons about 
their equilibrium position and their interaction with the moving fluxons. However the basic 
physics of the problem i.e. the shift of the voltage values is independent on the choice of the 
boundary conditions. 



III. LONG JUNCTION LIMIT 

A 4 th order Runge Kutta method with fixed time step At = 0.01, was used for the 
integration of the equations of motion. The number of grid points is N = 1000. We discuss 
first the case where the junction length is long Z = 20. We present in Fig. fj the I — V 
characteristics for the first and second ZFS that correspond to the case where one or two 
fluxons are moving into the junction. The pairing state is B\ g (a), E u (b), B\ g x E u (c). 
For the B\ g the external current cannot move the fractional ffuxon (//) which is confined 
at x = (see Fig. 0(b)). However for certain value of the bias current the // is transformed 
into an fractional antifluxon (faf) and an integer ffuxon (F) is emitted which is traveling 
to the left. The F hits the left boundary and transforms into an integer antifluxon (AF) 
which moves to the right. When the AF reaches the center interacts with the faf but is 
not able to change its polarity and results into an faf and an AF moving to the right. The 
antifluxon hits the right boundary transforms into a fluxon which moves to the center where 
it meets the oscillating faf and interacts with it forming a // and the period is completed. 

When a faf exists at the junction center, by applying the external current it emits an 
AF which moves to the right and it converts itself to a // (see Fig. 0(a)). The AF hits 
the right boundary and transforms into a F which moves to the left. When the F reaches 
the center interacts with the // and results into a // and a F moving to the left. The 
fluxon hits the left boundary transforms into an antifluxon which moves to the center where 
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it meets the oscillating // and interacts with it forming a faf and the period is completed. 

In the relativistic limit c$ ~ 1 reached at high currents a full period of motion back and 
forth takes time t = lljcs = 40 and since the overall phase advance is 4ir the normalized 
voltage will be V — <pt — 40 — 0.314 for the first ZFS. So when the junction length is large, 
the ZFS occur at the same values of the dc- voltage independently from the pairing symmetry, 
since one full fluxon or antifluxon propagates in the junction. These values much exactly the 
ones for conventional s-wave superconductors junction. The direction of the voltage pulse 
depends on the sign of the intrinsic flux and can be used for the qubit implementation. 

The different character of the various fluxon solutions can also be seen from the plot of the 
instantaneous voltage <pt at the center of the junction for the various fluxon configurations. 
This plot is seen in Fig. ^ for the solutions regarding the first ZFS. During the time of 
one period three peaks appear in this plot by the time when the fluxon (antifluxon) passes 
through the junction center. For the first ZFS the (j>t vs t plot can be used to probe the 
existence of // or faf at the junction center. The height of the middle peak is smaller 
for a bound // at the junction center, than for a bound faf. The plot of <pt at the edges 
shows two peaks at time instants which differ by half a period. Note that the characteristic 
oscillations of <p t between the peaks are due to the oscillation of the bound solution about 
the junction center. These oscillations have the same amplitude for the // and faf cases. 

For the Bi g case the bound fluxon and antifluxon have equal magnitudes or contain equal 
magnetic flux and the critical current is the same as seen in Fig. |2|(a). However for the E u 
case the // contains less flux than faf and has smaller critical current as seen in Fig. @(b). 
For the B ig x E u case the // contains more flux than faf and has greater critical current 
as seen in Fig. @(c). 

For the second ZFS multiple solutions exist in which two fluxons are propagating in 
the junction in different configurations. These solutions can be classified depending on the 
fluxon separation as seen in Fig. [| and give distinct critical currents in the I — V diagram. 
For all the modes in the second ZFS we can estimate the value of the constant dc-voltage 
where they occur as follows. A full period of motion back and forth takes time T = 40, and 
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since the overall phase advance is 8ir, in the relativistic limit where u = 1 reached at high 
currents, the dc voltage across the junction will be V = 0.628. So compared to the case of 
conventional s-wave superconductors junction we observe several curves for the second ZFS 
depending on the relative distance between the ffuxons and this may be used to probe the 
presence of intrinsic magnetic flux. 

We also considered damping effects due to the quasiparticles because the Josephson junc- 
tions made of high-Tc materials are highly damped. In the I — V curves higher damping 
shifts the I — V curve upwards (see Fig. [|(a)) and the fluxon reaches the critical current 
velocity only for currents that are very close to the critical current where the jump to the re- 
sistive branch occurs. In this case the ZFS are not 'vertical'. In the plot of the instantaneous 
voltage at the middle of the junction versus the time (Fig. 0(b)) the difference in height 
between the peaks when the F or the AF interacts with the bound // is larger for greater 
values of damping. The oscillations between the peaks become very large by increasing the 
damping and finally make the solution unstable. We believe that these solutions are more 
stable for small values of the damping at least for inline boundary conditions. In the overlap 
case (not presented in the figure) these solutions are more stable compared to the inline 
case. 

IV. SHORT LENGTH LIMIT 

When the junction length is small I = 2 the fractional fluxon or antifluxon does not 
remain confident at the junction's center but is able to move along the junction as a bunched 
type solution. The moving fluxon configuration could have fractional flux and additional 
steps are introduced in the I — V diagram. We plot in Fig. ^ the I — V characteristics for 
the Big, E u , Big x .E^-wave pairing states, for the small junction length 1 = 2. 

For the Bi g pairing state, and the first ZFS the moving fluxon configuration to the 
right is a combination of a fractional fluxon and a fractional antifluxon which contains half- 
integer magnetic flux (see Fig. 0(a)). When it hits the right boundary it transforms to 
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a configuration with opposite sign of flux which moves to the left. Then it hits the left 
boundary and transforms into a configuration that moves to the center and the period is 
completed. In the relativistic limit reached at high currents a full period of motion back 
and forth takes time t = lljcs = 4 and since the overall phase advance is 2n the normalized 
voltage will be V = ~ as seen in Fig. [](a) for the solution labeled as 1/2. Note that this 
value is half than the case where a full fluxon moves into the junction. In Fig. |g(a) we plot 
the <fit vs t at the center of the junction where the successive peaks correspond to the passage 
of the fluxon combination from the junction center. The (fit pulse is composed of two peaks, 
a positive and a negative one, corresponding to the AF and the // for the fluxon traveling 
in the forward direction. The pulse structures corresponding to the forward and backward 
directions are the same due to the symmetrical configuration of the fluxons traveling in the 
forward and backward directions. 

It is also possible to have solutions where an integer fluxon plus a fractional fluxon is 
propagating into the junction (see Fig. 0(b)). In this case the magnetic flux is equal to 
1.5. In a junction of length I = 2 the propagating fluxon accomplishes an overall phase 
advance of 67r in a full period T = 4. Thus the voltage across the junction will be V = 37r/2 
as seen in Fig. |B|(a) for the solution labeled as 3/2. In Fig. ||(b) we plot the (fi t vs t 
at the center of the junction where the successive peaks correspond to the passage of the 
fluxon combination from the junction center. The moving fluxon or antifluxon has internal 
structure and therefore a double peak structure appears in the (fi t vs t diagram. 

Finally in Fig. 0(c) we present 0(x), for the case where two integer fluxons with a 
fractional fluxon in between move into the junction. This configuration contains magnetic 
flux equal to 2.5. Thus the dc voltage across the junction will be V = 5ir/2 as seen in Fig. 
^|(a) for the solution labeled as 5/2. In Fig. [|(c) we plot the <fi t vs t at the center of the 
junction where the periodic pattern of three peaks corresponds to the passage of the integer 
fluxons and the half integer fluxon from the junction's center. 

So for the B\ g pairing state the ZFS appear at values of the normalized voltage that 
are displaced by 1/2 compared to the case of s-wave superconductor junction. This value 
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of the voltage corresponds to the intrinsic phase shift. This is analogous to the shift of the 
critical current versus the magnetic flux in a corner squid of ci-wave superconductors || and 
is expected to be confirmed by experiment. 

For the E u and B\ g x E u cases, we can have additional solutions where the moving fluxon 
has integer flux and the voltage steps appear at values V = nn in addition to V = ^ as 
seen in Figs, ^(b) and |6](c) for the solutions labeled with integer numbers. For the _Bi 9 -wave 
case the forward and backward configurations are symmetric and successive peak structures 
have the same form. In the E u and B\ g x E u cases, successive peaks, corresponding to the 
structure of the fluxon configuration that moves in the forward and backward direction, 
have different amplitudes indicating that the fluxon configurations moving in the forward 
and backward directions have different structure. 

We examined also the case where the pairing symmetry of the superconductor is d + is. 
In this case due to the difference in the flux content of the static solutions [|J, the critical 
currents for the //, faf modes of the first ZFS do not coincide. 



V. QUBIT IMPLEMENTATION 

Also the frustrated junction could be considered as a way to build a qubit. This idea 
has also been implemented using s-wave/ci-wave/s-wave junction exhibiting a degenerate 



ground state and a double-periodic current-phase relation |21| , or superconducting Josephson 
junction arrays |22fl . Also the dynamics of a Josephson charge qubit, coupled capacitively 
to a current biased Josephson junction has been studied f23fl . The two segments of the 
frustrated junction have characteristic energies E(<p), and E(<j) + 7r) and the resulting energy 
has minima at ±7r/2. The system exhibits a degenerate ground state. The bound //, faf 
can be considered as the two quantum levels of our system. By applying an external current 
a bunched type solution containing 0.5(— 0.5) flux is generated if the actual ground state of 
the system is ff{faf) which propagates into the junction to the left(right) and generates 
a voltage pulse which can be determined at the boundaries. A possible method to generate 
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the desired ground state for our system, i.e., //, or faf is to apply an external magnetic 
field (positive or negative) and then to slowly decrease it. Depending on the sign of the 
external field the system will go either in the // or in the faf state. 

In the E u , B\ g x E u cases the actual ground state of the system is not degenerate, i.e., the 
// and faf carry different flux. Moreover the traveling fluxon caring half the flux quantum 
to the left has different structure than the one traveling to the right and the <f>t at the ends 
are different for the //, faf. In the long length limit the presence of the //, faf at the 
center can be deduced from a measurement of the <f>t at the edges since the direction of the 
moving integer flux depends on the sign of the intrinsic flux. In this case the intrinsic flux 
can not escape from the junction's center. 



VI. CONCLUSIONS 

The fluxon dynamics in frustrated Josephson junction with p, d, and /-wave pairing 
symmetry, is different in the long and short junction limit. When I is large, the bound 
intrinsic flux remains confined at x — 0, and the moving integer fluxon or antifluxon interacts 
with it only when it approaches the center. However when the length is small the bound 
fluxon becomes able to move as a bunched type solution. For <i-wave junction the I — V 
curves are displaced by a voltage that corresponds to the intrinsic phase shift. 

The resonant fluxon motion also can be determined experimentally in one-dimensional 
ferromagnetic — 7r junctions where the width of the ferromagnetic oxide layer determines 



the region of the junction where the Josephson critical current is positive or negative p4| , |25 
In this case the junction contour does not have to change as in the case of a corner junction. 
However the change of the junction contour was not a problem for the realization of the 
intrinsic fluxon in the static problem, so way should it prevent the realization of its motion? 

A final comment is that the frustrated junctions that we consider in this paper are 
realized in the (a6)-plane of the unconventional superconductors due to the sign change of 
the order parameter. This type of junctions is different from the series array of intrinsic 
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Josephson junctions in high-Tc supercondcutors where the Josephson effect is observed in 
the c-axis, for instance in Bi 2 Sr 2 CaCu20 8 crystals p6| . 
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(a) 
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a 



(b) 



d,p, or f-wave 




1 




s-wave 



FIG. 1. (a) A single Josephson junction between superconductors A and B with a two com- 
ponent order parameter. Also a small co-ordinate system indicating a and b crystalline directions 
is shown, (b) The geometry of the corner junction between a mixed symmetry superconductor and 
an s-wave superconductor. 
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FIG. 2. Normalized current versus the normalized voltage for the inline geometry, for the first 
and second ZFS. The solutions for the first ZFS are the //, faf corresponding to a bound fluxon 
or antifluxon in the junction's center. For the second ZFS the solutions are labeled by the relative 
distance between the fluxons l/x, where I = 20 is the junction length, and x = 1, 2, 3, 6 respectively, 
7 = 0.01. The pairing state is (a) i?i 5 -wave, (b) E u -w&ve, (c) B\ g x E u -w&ve. 
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FIG. 3. Phase 4>(x) vs x for the solutions in the first ZFS, at various instants, during one 
period separated by At = 2.8. The curves are shifted by 0.5 to avoid overlapping. I = 20, I = 1.6, 
7 = 0.01: (a) //, (b) faf. The pairing state is i?i s -wave. 
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FIG. 4. Instantaneous voltage in the middle of the junction (x = 0) vs time t, for the solutions 



in the first ZFS. / = 20, 7 = 0.01, / = 1.6: (a) //, (b) faf. The pairing state is B\ g -wme. 
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FIG. 5. (a) Normalized current versus the normalized voltage for the inline geometry, for the 
// solution, for different values of the damping constant 7. I = 20. The pairing state is Bi g -wave. 
(b) Instantaneous voltage in the middle of the junction (x = 0) vs time t, for the solution //, 
/ = 20, / = 1.6. The pairing state is £>i g -wave. The damping constant is 7 = 0.03. 
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FIG. 6. Normalized current versus the normalized voltage for the inline geometry, for a junction 
of length I = 2, 7 = 0.01. The different modes are labeled by the value of the normalized voltage 
divided by n in which they occur. The pairing state is (a) i?i 9 -wave, (b) E u -wave, (c) B\ g x E^-wave. 
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FIG. 7. Phase <j)(x) vs x for the solutions in the first three ZFS, at various instants, during one 
period separated by At = 0.2. The pairing state is £>i 9 -wave. The curves are shifted by 0.5 to 
avoid overlapping. I = 2, I = 0.25, 7 = 0.01. 
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FIG. 8. Instantaneous voltage in the middle of the junction (x = 0) vs time t, for the solutions 
in the first three ZFS. The pairing state is Bi g -wave. I = 2, 7 = 0.01, / = 0.25. 
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